Importance of matrix elements in the ARPES spectra of BISCO 
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We have carried out extensive first-principles angle- resolved photointensity (ARPES) simulations 
in Bi2212 wherein the photoemission process is modelled realistically by taking into account the 
full crystal wavefunctions of the initial and final states in the presence of the surface. The spectral 
weight of the ARPES feature associated with the Cu02 plane bands is found to undergo large and 
systematic variations with fc|| as well as the energy and polarization of the incident photons. These 
theoretical predictions are in good accord with the corresponding measurements, indicating that 
the remarkable observed changes in the spectral weights in Bi2212 are essentially a matrix element 
effect and that the importance of matrix elements should be kept in mind in analyzing the ARPES 
spectra in the high-Tc's. 
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Angle-resolved photoemission spectroscopy (ARPES) j/^ i^^\ 
has contributed significantly towards an understanding 
of the nature of the normal as well as the superconduct- 
ing state of the cuprates |l]-0] . Much of the data on the 
cuprates, however, has been analyzed by assuming that 
the ARPES essentially measures the one-particle spectral 
function of the initial states. While this simple approach 
yields insights into the underlying physics, a satisfactory 
description of the spectra must necessarily model the 
photo-excitation process properly by taking into account 
the matrix element involved, the complex modifications 
of the wavefunctions resulting from a specific surface ter- 
mination, and the effects of multiple scattering and of 
finite lifetimes of the initial and final states. 

This article presents first-principles simulations of 
ARPES spectra in Bi2212 (BISCO) using the one-step 
model of photoemission which incorporates the afore- 
mentioned effects realistically p-M. We focus on the 
ARPES signature of Cu02 plane bands which are widely 
believed to be the key to the mechanism of supercon- 
ductivity in the cuprates. The spectral weight of the 
ARPES peak associated with the Cu02 plane bands is 
found to undergo large variations with kn as well as the 
energy and polarization of the incident photons. These 
theoretical predictions are in remarkable accord with the 
corresponding measurements |11[| on BISCO, and show 
clearly the importance of " matrix element effects" |1^ in 
the ARPES spectra. Notably, a substantial increase in 
the ARPES spectral weight in BISCO in going from F 
to M was noted early by Anderson |1^], who speculated 
that this puzzling behavior may be the hallmark of spin 
charge separation. 

The physical and formal underpinnings of our ap- 
proach may be exposed by starting with the following 
Golden rule-based expression for photointensity from ini- 
tial states at energy E with photons of energy Tiuj: |Q P{,E, lo) 



^i(^y) are the initial (final) states of the semi-infinite 
solid, and A = e'h/2mc{p • A -I- A • p) is the interac- 
tion Hamiltonian with the electron momentum operator 
p and the vector potential A of the photon field. 

In the one-step model used in the present computa- 
tions, Eq. 1 is manipulated into the form H 

im,E,hu;) = --Im < ku\G^ {E + hLj)AGt iE)A'' 
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where the matrix element involves the free electron final 
state of momentum kn . G2 and Gi denote the retarded 
(+) or advanced (-) one-electron Green functions at ap- 
propriate energies. 

Notably, the so-called three-step model of photoemis- 
sion approximates the matrix element in Eq. 1 via the 
bulk Bloch wavefunctions yielding for the photointensity 
within the solid 115(1 : 
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which is cast in terms of the one-particle spectral func- 
tions of the initial and final states {Ai and Af), and the 
processes of transport and emission are to be considered 
separately. Assuming that: (i) The system is strictly 2D, 
and (ii) the final states form a structureless continuum, 
Eq. 3 reduces to: 

~ " (4) 
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in terms of just the Ai{E)'s. Further, for a single band 
solid, one obtains: 
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where the Fermi function on the right side is suppressed. 

Brief comments on forms 1-5 are appropriate in order 
to highUght the underlying approximations. Since forms 
3-5 ignore the presence of the surface, it is very difhcuh to 
include effects of different surface terminations n& , and 
the associated distortions of bulk wavefunction which can 
be quite severe. In form 4, even near an ARPES peak 
from a specific initial state, other states will in general 
contribute a background upon being broadened due to 
their finite lifetimes. Forms 4 and 5 ignore the struc- 
ture in the final state spectrum which will be seen below 
to be quite significant. We emphasize that the distinc- 
tions between the processes of excitation, transport and 
emission through the surface invoked by the three-step 
model (Eqs. 3-5) are artificial since the more satisfac- 
tory one-step formula (Eq. 2) does not admit such a 
decomposition. 

The relevant technical details of our computations are 
as follows. In order to keep the problem manageable, 
the modulation of the lattice is neglected and the crystal 
structure of BISCO is assumed to be perfectly tetragonal; 
this still involves 30 atoms per unit cell, and a substan- 
tial extension of the earlier work on relatively simpler 
lattices 0,0 The crystal potential was first obtained 
self-consistently within the KKR scheme [|l^ , and essen- 
tially yielded the well-known LDA-based band structure 
and Fermi surface of Bi2212 |l8]; however, the actual 
potential used here is slightly modified in that the Bi-0 
pockets around the Af -point are lifted above the Fermi 
energy to account for their absence in the experimen- 
tal spectra. |13] The surface is assumed to terminate in 
the Bi-0 layer in accord with the general consensus pH] 
that this is an easy cleavage plane in BISCO; the inci- 
dent light is assumed polarized in the {k^, ky) plane since 
the ARPES spectra of interest here seem to be insensi- 
tive to the z-component A^ of the vector potential, and 
further, dielectric effects would in general complicate the 
relationship between the value of A^ inside and outside 
the crystal Q. 

The ARPES simulations of Fig. 1 help set the stage for 
our discussion. As expected, the separation between the 
spectral features A and B in Fig. la associated with the 
Cu02 plane bands increases as we go from F towards M . 
p2| Figs, lb and Ic show that the spectral weight of A or 
B depends dramatically on the values of kn and photon 
energy hu. The weights display large changes for either 
peak as a function of hv for a fixed fcy , or as a function of 
kn for a fixed hv. These results highlight the importance 
of matrix element effects in BISCO since the weights will 
be constant independent of hv or kn in the absence of 
these effects. In fact, our theoretically predicted hv and 
kn dependencies of spectral weights are in substantial ac- 
cord with the measurements; we discuss this aspect now 
with reference to Figs. 2-4. 

Fig. 2 compares the measured and computed total 
spectral weights in a 500 meV energy window below Ef 



along the three high symmetry directions. The charac- 
teristic features of the experimental data are: a relatively 
low flat intensity along TX, a steep rise along TM com- 
pared to ry, and the presence of the prominent bump 
around A:|| ~ 0.25 along TY . All these aspects of the data 
are essentially reproduced by the theory. ||2^ 

Fig. 3 considers the photointensity in the {kx,ky) 
plane for two different polarizations of the incident light 
where the initial state is held fixed at Ef. The color 
plots give computed intensities over a dense grid of fcy- 
values and display the two Cu02 plane band sheets A 
and B in the band structure; these plots are representa- 
tive of what could be measured in a suitably arranged 
constant-initial-energy angle-scanned (CIE-AS) ARPES 
measurement. gj2j,^ Incidently, Figs. 3a2-3c2 are not 
symmetric about a horizontal or vertical line through the 
center. In Figs. 3b2 and 3c2, the light is incident along 
the TY direction and the figure therefore is symmetric 
only about this diagonal line. In Fig. 3a2, on the other 
hand, the light is polarized horizontally, and the inten- 
sities are symmetric around the TM line; this symmetry 
becomes visible only when the figure is extended in the 
vertical direction to include a larger range of momenta. 

The intensity associated with the outer plane band 
sheet A as one goes around the Fermi surface is shown in 
quantitative detail in Figs. 3al-3a3. Since A is generally 
more intense than the inner plane band B, A is presum- 
ably more relevant in connection with the experimental 
data near Ef. For light polarized along the horizontal 
direction (Fig. 3a), we see that the intensity is large 
around a = 0° and decreases rapidly beyond a k, 20°. 
A 45° rotation of the polarization vector (Figs. 3b and 
3c) induces substantial changes in the shape and magni- 
tude of the intensity, and the appearance of a minimum 
at /? or 7 = 45°. The experimental points in Fig. 3 are 
in good overall accord with the measurements, some dis- 
crepancies around a and /3 ~ 0° notwithstanding. Thus 
quite large observed variations (nearly an order of mag- 
nitude) in the emission intensity from different parts of 
the Fermi surface are mainly the consequence of matrix 
element effects. 

Finally, Fig. 4 considers the photon energy dependence 
of the emission intensity of the spectral feature around 
the M-point. The theoretical curve displays a prominent 
peak around 22 eV and indicates clearly that the final 
states in BISCO possess considerable structure which is 
neglected in approximations of Eqs. 5 and 6, especially 
when the matrix element in these equations is further 
replaced by a constant. The experimental points which 
show the presence of a broad peak centered around 21 
eV are in good overall agreement when we keep in mind 
that errors of the order of a few eV's in locating the final 
states are generally expected in the first principles band 
structures. 

We emphasize that the agreement seen in Figs. 2-4 is 
robust to uncertainties inherent in such a comparison on 



the theoretical as well as the experimental side. The ex- 
perimental weights depend on the specific energy window 
used in their definition, and on whether or not a suitable 
background is subtracted. |£6[ In this vein, the computed 
weights depend also of course on the specific values of the 
initial and final state damping parameters E^ and S ^ . In 
order to assess these effects, we have carried out exten- 
sive simulations using a variety of different values and 
energy-dependencies of S^ and S p in addition to vary- 
ing the real part of the initial state self-energy (in order 
to mimic correlation effects, even though the LDA frame- 
work underlies our computations) , and find that the main 
features of the results of Figs. 2-4 are insensitive to such 
variations. 

In summary, we have carried out extensive first- 
principles one-step ARPES simulations in Bi2212 
wherein the photocmission process is modelled realisti- 
cally by taking into account the nature of the initial and 
final state crystal wavefunctions as well as the multi- 
ple scattering effects in the presence of a specific sur- 
face termination. We focus on the nature of the ARPES 
feature arising from CUO2 plane bands, and consider 
in particular its spectral weight as a function of fcy as 
well as the energy and polarization of the incident pho- 
tons. Large variations in the spectral weights predicted 
theoretically along three different high symmetry direc- 
tions are in good accord with the corresponding mea- 
surements. A good agreement between theory and ex- 
periment is also seen with regard to changes in spectral 
weights with photon energy around the M-point, as well 
as along the Fermi surface contours in the {kx,ky) plane 
for two different polarizations. This study shows clearly 
that the remarkable observed changes in the ARPES 
spectral weights in Bi2212 are essentially a matrix ele- 
ment effect and that the importance of matrix elements 
should be kept in mind in analyzing the ARPES spec- 
tra in the high-Tc's. Another notable implication of this 
work is that the integral (over energy) of the ARPES 
intensity does not yield the momentum density of the 
electron gas. 
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Figure captions 

Figure 1: (a): Simulated ARPES intensity in BISCO 
for k\\ varying between T and A/ (bottom to top) at 
hv = 22ey. (b) and (c) show computed variations in 
spectral weights (areas under peaks) of the two Cu02 
plane band features A and B with photon energy for dif- 
ferent /c|| -values, based on a crystal potential without Bi- 
O hole pockets at M. A small value of the initial and 
final state damping parameters is used to highlight spec- 
tral features. 

Figure 2: Theoretical weights obtained by integrating 
the hv — 22ey ARPES spectra over the binding energy 
range of 0-500 meV are compared with the corresponding 
experimental results |0]. Curves for the three symme- 
try lines (Brillouin zone in inset) are offset vertically for 
clarity; kn is given in relative units such that the distance 
from r to X, Y or M is defined to be unity for each di- 
rection. Theoretical values are normalized to match the 
experimental value at the maximum around A:|| w 0.8 in 
the TM curve. 

Figure 3: The color plots give simulated ARPES in- 
tensities for emission from Ef at hv — 22 eV for two dif- 
ferent polarizations (white arrows) of the incident light. 
The Cu02 plane band sheets are marked A and B. The 
intensity (area under peak) of the outer plane band A is 
shown in (al)-(cl) as a function of the angles a, /3 and 
7. Experimental data after Ref. ||T^. Theory normalized 
around a w as shown in (al). 

Figure 4: Spectral weight (over a 500 meV window) of 
the feature at Af -point is compared with the correspond- 
ing experimental ||ll[ results as a function of the photon 
energy. Theory normalized to experiment around 21 eV. 
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